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Abstract 

The exclusive deuteron break-up reaction is analyzed within a covariant ap- 
proach based on the Bethe-Salpeter equation with realistic meson-exchange inter- 
action. Relativistic effects in the cross section, tensor analyzing power and polar- 
ization transfer are investigated in explicit form. Results of numerical calculations 
are presented for kinematical conditions in forthcoming p + D reactions at COSY. 
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1. Introduction: Break-up reactions of deuterons by protons receive presently a renewed 
interest |TJ, 0, [3], M . The investigations of break-up processes are motivated by the hope 
to extract directly the deuteron wave function from experimental data, supposed the 
mechanism described by the impulse approximation dominates at moderate values of 
intrinsic momenta. Pioneering experimental studies of elastic backward and inclusive D 
+ p reactions have been performed in Dubna and Saclay || H |J . It is found, indeed, that 
the impulse approximation holds in a large interval of momenta of detected protons, except 
for a broad shoulder at momentum of the outgoing proton (measured in the deuteron's 
rest frame) of about 0.3 GeV/c. Probably small corrections to the impulse approximation 
are sufficient to account for this shoulder. In the same experiments the polarization 
observables of the deuteron, such as the tensor analyzing power T20 and polarization 
transfer k are measured. These quantities are more sensitive to the reaction mechanism, 
and a combined analysis of them may provide information about S and D components 
of the deuteron wave function separately. However, as in the unpolarized case, the data 
on T20 and k exhibit systematic deviations from theoretical predictions |5|, ||] in the same 
momentum region of the outgoing proton. A kinematical analysis of the invariant missing 
mass shows that this region corresponds to A excitations, so that additional corrections 
should supplement the calculations based on the impulse approximation. From this it 
becomes clear that an interpretation of the inclusive data in terms of a deuteron one- 
body momentum density becomes ambiguous. Furthermore, the typical energies in these 
reactions seem already too high for a non-relativistic approach. It is necessary, therefore, 
to describe these processes in a covariant formalism. 

The measurement of proton-deuteron break-up reactions with polarized particles in an 
exclusive experimental setup is planned at the cooler synchrotron COSY in Jiilich [|IJ. 
The experiment will detect the scattered fast proton in the forward direction in the spec- 
trometer ANKE in coincidence with the slow backward-emitted proton. The kinematical 
conditions are chosen in such a way that the missing mass of the undetected particles 
is exactly equal to the neutron mass. Hence, this experiment offers a unique possibility 
of using polarized deuterium targets in combination with a polarized proton beam for 
extended studies of exclusive deuteron break-up processes. The coincidence measurement 
allows one to exclude kinematically particle production processes and, consequently, to 
constrain the reaction mechanism. Two important aspects in these experiments ought to 
be stressed here: (i) The possibility of excluding A excitation processes will allow for un- 
derstanding better the role of the meson degrees of freedom and to clarify the origin of the 
above mentioned shoulder in the inclusive processes, (ii) A comparison with theoretical 



predictions will provide a good test of the validity of the spectator mechanism. 

Previous relativistically invariant investigations |S], [7]] of one-nucleon exchange diagrams 
for inclusive and elastic D + p scattering are based on numeric solutions of the Bethe- 
Salpeter (BS) equation with a realistic interaction. Several authors || studied relativistic 
effects in the deuteron by considering the D — > NN vertex within approximations to the 
exact BS equation and analyzed so electromagnetic and elastic hadron scattering off the 
deuteron. Up to now, however, a consistent investigation of polarization phenomena in 
exclusive processes with explicit identification of relativistic effects is still lacking. 

In the present work we perform such a covariant analysis of the exclusive proton - 
deuteron break-up reaction. Fully covariant expressions for the cross section, tensor an- 
alyzing power T20 and transferred polarization k are obtained within the BS formalism. 
The contribution of pure relativistic corrections is separated. Results of numerical cal- 
culations, utilizing the recently obtained numerical solution of the BS equation with a 
realistic interaction kernel, are presented for kinematical conditions available at COSY. 
2. The spectator mechanism: We consider the above described exclusive break-up 
reaction of the type p + D — > pi(0°) + p2(180°) + n3(0°). Within the spectator mechanism 
approach this reaction is presented by the Feynman diagram shown in fig. la, where the 
upper and lower vertices factorize and, consequently, they can be computed separately. 
According to the above kinematical conditions one may consider the fast forward proton 
to be produced by elastic scattering of the beam proton off one nucleon in the deuteron. 
In fig. lb this part of the diagram is depicted and the corresponding cross section reads 
,gnn 1 5 4 (p p + p n -p 1 -p 3 ) 2 d 3 Pl d 3 p 3 

dB W 5 W 

where X(p p ,p n ) = (PpPn) 2 — m 1 r n n an d the invariant nucleon-nucleon (NN) amplitude Tnn 
is expressed via the truncated vertex functions T NN and projectors A aa > as follows 

\Tnn\ 2 = - A-aa'iPp) A/3/3>(p n ) O a j3]a' f3' ■, 

O a p- a if3i = ( M G°3, S 3 )T NN u(p U Si)) (u(pi, S!)r NN u(p3, S 3 )) a/3 , (2) 

spins 

where u denote nucleon Dirac spinors and a, (3 are Dirac indices (with summation over 
twofold indices) and the operator @ is related to the T matrix of nucleon-nucleon scat- 
tering. To calculate the latter one in a fully covariant manner one needs an analysis of the 
NN vertices within an appropriate relativistic formalism, e.g., the BS approach. For the 
elastic scattering the NN amplitude is thoroughly investigated in ref. ||. In our case the 
operator O determines the differential elastic cross section da/dt at low transferred mo- 
mentum t = (p p — pi) 2 , and the contribution of the spin part in eq. (§) may be neglected. 



Then the operator O a p- a rpi takes the simplest possible form O a p. ia ip> ~ 8 aa '5f3f3/\A NN \ 2 , 
where Ann is the scalar part of the NN amplitude. Then, with m as nucleon mass, 

da NN = -—^——\A NN \ 2 dt. (3) 

Now we are in the position to calculate the main contribution to the process described by 
the diagram in fig. la. Applying the Mandelstam method the invariant cross section for 
the break-up of a polarized deuteron with spin projection M can be written as 



,2 m 2\ 



r/'rr 1 jMjWPn) da NN (pj -///-;.,, T n , irD .,.. ,\ 

dtd 3 p 2 (2tt) 3 ^\(p p ,p D ) dt 2m V 7 

where p = (p n — p2)/2, and / stands either for a unity matrix in case when the outgoing 
proton p 2 is not polarized or for (1 + 75S2V2 otherwise (s 2 is the contracted polarization 
four vector). ^ffj(p) denotes the charge conjugated BS amplitudes introduced in ref. |I0|; 



numerical solutions are reported in refs. [11, 12 



Eq. (f|) is a rather formal result and an interpretation of different contributions is 
straightened in the present form. A few subtle aspects are to be mentioned: 

(i) The detected proton is on mass-shell so that p\ — m? = 0, and the cross section seems 
to vanish. However, the BS amplitude itself is singular when one nucleon is on mass-shell. 
To get a finite result the corresponding expressions should be evaluated analytically. 

(ii) The numerical solution of the BS equation is obtained in the Euclidean space-time, 
where the time component po of the relative momentum p is purely imaginary. In processes 
under consideration p is fixed and real. Hence, one needs either a numerical procedure 
for an analytical continuation of the amplitudes to the real relative energy axis (cf. |7]]) 
or another recipe |13| for using our numerical solutions for this case . 



(iii) In solving the BS equation we expand the amplitude on the complete set of the 
Dirac matrices and obtain eight partial amplitudes for ]IU| , [TJJ]. However, it is known 
J7L that in such cases, where one nucleon is on mass-shell, only four partial amplitudes 
contribute to the deuteron observables. From eq. (|4]) it is not clear which amplitudes play 
the most important role in the process. 

To tackle the above problems it is convenient to transform our representation of the 
partial amplitudes to the so-called p spin classification |TJ|, [19. In this representation 
one usually adopts the spectroscopic notation for the partial amplitudes ( 2S+1 ) Lj lp2 . (For 
the explicit form of the unitary transformation matrix between these two representations 
cf. ||14|| .) In this notation it becomes immediately clear that in processes with one nucleon 
(say the second one) on mass-shell the relevant contribution to the cross section comes 



only from four amplitudes with positive second p spin index, i.e., 3 Sf + , 3 Df + , l P{ + and 
3 P{ + for which we use the short hand notation \1/ 5, typ, '$p 1 and \l/p 3 , respectively. 
3. Relativistic corrections: The trace in eq. is evaluated by an algebraic formula 
manipulation code which delivers the contribution to the unpolarized cross section 
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where another zero due to 2E 2 +2p — Mp, = appears explicitly in the cross section (Mp, is 
the deuteron mass). To handle these two zeros and the singularities in the amplitudes 
it is convenient to introduce instead of the corresponding partial vertices G( P o, |p|) 
that have no poles when one particle is on mass shell. For an explicit relation between 



partial amplitudes and partial vertices we refer the interested reader to ref. [14]], where 
the dependence of S and D wave vertices upon the relative energy is shown to be smooth, 
contrary to the amplitudes which display a strong dependence on p . Therefore, in our 
calculations, we can replace, at moderate values of Pq, the S and D vertices by their 
values at po = with good accuracy. The P vertices can be expand into Taylor series 
about po = up to a desired order in po/m. Then the corresponding derivatives can 
be computed numerically along the imaginary axis since they are analytical functions of 
Pq |[HJ. After replacing the amplitudes by vertices, the zeros and singularities cancel at 
P20 — > E 2 , and the result is finite. Finally, to cast our formulae more familiar form, known 
from non-relativistic calculations, we introduce the notion of BS wave functions [0, ||, [14 



u(\p\) = f 



Gg(0,|p|) 
2E 2 - M D ' 



w n h _ Gg(0,|p|) 

Wm ~ :F 2E 2 -M^ Vp ^ 



(IpI) 



(o, |p|) 



M. 



D 



(6) 



where T = l/4iiy/2M^. With these definit ions one gets 
(i) the unpolarized differential cross section: 

M D \]KPp,Pu) da NN 
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(7) 
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(ii) the tensor analyzing power T 2 q: 

-W 2 (\p\) - 2v / 2t/(| P |)^(|p|) + 2V£(|p|) - V| 3 (|p|)l (9) 



IY, D m(\p\)) V2T : 
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(10) 



(iii) the polarization transfer k: 

(^E^(IpI))* = 



t/ 2 (lpl)-^ 2 (lpl) + ^f/(lpl)^(lpl) 
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(11) 
(12) 



For short hand notation we introduce the deuteron structure factor D^(|p|) with a defi- 
nition which follows from eqs. ([7[ §). 

A numerical analysis of solutions of the BS equation in terms of amplitudes within the 
p spin basis shows |nj that the BS wave functions t/(|p|) and iy(|p|) are the dominant 
ones and coincide to a large extent with the corresponding non-relativistic wave functions 
found as solutions of the Schrodinger equation with one-boson-exchange potential. The 
remaining two functions Vp 13 (|p|) are a few orders of magnitude smaller, and for the 
considered processes with momentum of the backward proton | p»2 1 = |p| < 0.51 GeV/c 
the diagonal terms in Vp 13 (|p|) are negligible. Therefore, eqs. (|7|, |9], [11]) are identified as 
the main contributions to the corresponding observables and they might be compared with 
their non-relativistic analogues. The interference terms (§, |10|, 0) posses contributions 
from negative states and are proportional to |p|/m. Due to their pure relativistic origin 
we refer to them as relativistic corrections in the deuteron break-up reactions. Note that, 
when disregarding the relativistic corrections and equating the wave functions Z7"(|p|) 
and VF(|p]) to their non-relativistic analogues, our formulae (0 - [12]) exactly recover the 
non-relativistic expressions for a, T 20 , k computed within the spectator mechanism ||. 
4. Results and discussions: In figs. 2 - 4 we present the results of our calculations 
by exploiting the numerical solutions [O, [12|, |14| of the BS equation with a realistic 



interaction kernel with ir, u, p, a, rj, 5 exchange (parameters as in table 1 in ref. [JTXJ] ) . The 
unpolarized cross section computed by eqs. (0, ||) is displayed in fig. 2, where the flux 
factor has been put to unity. This flux factor reflects only the dependence of the cross 
section on the beam energy within the spectator mechanism approach. By disregarding 
it one obtains the energy independent part of the cross section. Within the spectator 
mechanism, both T 2u and k do not depend on the initial energy. For the elastic neutron 
- proton scattering we use a fit of data 0]. The dashed curves in figs. 2-4 depict the 



contribution of only positive waves according to eqs. (|7], |9|, |TT]), while the short-dashed 
curves are the relativistic corrections according to eqs. |X0, |TJ|) (in fig. 2 we display the 
modulus of the corrections, because of a sign change), and the full lines are their sums. 
For completeness we also present the results of non-relativist ic calculations with the Bonn 
potential wave functions shown as dotted curves. 

It is seen in figs. 2-4 that the relativistic corrections are negligible small in a wide 
range of momenta | p>2 1 and become significant only at | p>2 1 > 0.6 GeV/c. From this we 
conclude that at kinematical conditions as envisaged in COSY experiments Jl| the rela- 
tivistic corrections may be safely neglected since they are much smaller than the expected 
experimental errors. At E p ~ 3 GeV one has |p2|maa; ~ 0.51 GeV/c. (Note that within 
the spectator mechanism the maximum value of | j>2 1 is restricted to |p2|maa; ~ 0.8 GeV/c.) 
Hence, in the proposed experiments one may investigate in great detail different aspects 
of the reaction mechanism, or the contribution of non-nucleonic degrees of freedom like 
meson-exchange currents and the role of A isobars, without bothering about relativistic 
effects. Note that our expressions for the observables in exclusive break-up processes are 
similar to those for the inclusive ones |5], ||. There exist an essential advantage in the 
proposed exclusive experiments: in inclusive processes the detected shoulder in the 
cross section at | p»2 1 ~ 0.3 GeV/c cannot be explained within the spectator mechanism. 
The origin of the discrepancy is believed || to root in the contribution of meson produc- 
tion in the NN vertex. In the exclusive processes these contribution may be separated 
kinematically and conclusions about this problem can be settled. 

The relativistic corrections eqs. (|| [H| [12]) are governed by negative P states in the 
deuteron. This can be considered as a hint that admixtures of P waves within the BS 
approach are related to relativistic corrections by taking into account meson-exchange 
currents and NN pair production diagrams |T7]] in the non-relativistic picture. To establish 
a correspondence between our results and the mentioned non-relativistic calculations we 
estimate the contribution of the relativistic corrections by computing the P wave vertices 
in the so-called "one-iteration approximation" . The gist of this approximation is as follows 
[|T8]]: in solving the BS equation by an iteration procedure one puts as zeroth iteration 
the exact solution of the Schrodinger equation for S and D vertices and zero for other 
waves; then the P vertices are found by one iteration of the BS equation. Our experience 
in solving numerically the BS equation shows that it converges rapidly for relatively small 
momenta < 1 GeV/c. That means when utilizing the exact non-relativistic solutions, 
after one iteration the resulting P waves are not too far from the reality. 

Skipping cumbersome algebraic manipulations the result for the function V in eq. (^ 



with a BS kernel with pseudo-scalar one-boson exchange reads as follows 

Vp 1i3 (|p|) = -^ 2 tjV dr (1 + A*r) ji(|p|r) [N u U(r) + N w W(r)} , (13) 

where U(r) and W(r) are the non-relativistic deuteron wave functions in the coordinate 
representation, and g\ ~ 14.5 is the pion-nucleon coupling constant; N u = 1 (y/2) and 
N w = y/2 (— 1) for Pi (P3) waves. Introducing this result into the expression for the cross 
section eq. (0) the relativistic corrections is 

^ dV- c - \]KVp,Vn) da NN ( 8 v gl f J e~» r H . . . . 
E 2 - = = v , — ; — \—-— dr (l + ur)h(\p\r) 



x {[/(|p|)[-[/(r) + 2y/2W(r)] + W(\p\)[2y/2U(r) + W{r)}} (14) 

which is similar to expressions obtained in non-relativistic evaluations of the so-called 
"catastrophic" and pair production diagrams in electro-disintegration of the deuteron 
[pTgfl . In eq. flU]) v/c is the velocity of the detected slow proton, and the quantity in the 
large parenthesis may be interpreted as effective number of NN pairs in the deuteron 
contributing to break-up reactions (details will be presented elsewhere [0). From this 
it becomes clear that generic relativistic calculations, even in impulse approximation, 
contain already some specific meson-exchange diagrams, i.e., pair production currents, and 
one should pay attention on the problem of double counting when computing relativistic 
corrections beyond the spectator mechanism. 

5. Summary: In summary, we present for the first time an explicit analysis of relativistic 
effects in exclusive deuteron break-up reactions within the Bethe-Salpeter formalism with 
realistic interaction kernel. Numerical estimates of relativistic effects in the cross section, 
tensor analyzing power and polarization transfer at kinematical conditions of forthcoming 
COSY experiments are performed. Relativistic corrections under these conditions are 
identified and found negligible. It is shown that the planned experiments can discover 
mainly effects related to processes beyond the impulse approximation. 
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Figure captions 



Fig. 1: Feynman graphs for the exclusive reaction p + D = pi(0°) + p2(180°) + ^(0°) 
(a) and for the elementary NN vertex (b). 

Fig. 2: The spin averaged differential cross section E da/dtd 3 p 2 (with flux factor put to 
1) for the exclusive proton - deuteron break-up reaction. The meaning of the curves are 
explained in the text. 

Fig. 3: The deuteron tensor analyzing power T 2 q for the exclusive proton - deuteron break- 
up reaction. Notation as in fig. 2. 

Fig. 4: The polarization transfer k for the exclusive proton - deuteron break-up reaction. 
Notation as in fig. 2. 
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Fig. 1: Feynman graphs for the exclusive reaction p + D = pi(0°) + p2(180°) + n 3 (0°) 
(a) and for the elementary NN vertex (b). 



1 1 




Fig. 2: The spin averaged differential cross section E da / dt d 3 p 2 (with flux factor put 
to 1) for the exclusive proton - deuteron break-up reaction. The meaning of the curves 
are explained in the text. 
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Fig. 3: The deuteron tensor analyzing power T 2 o for the exclusive proton - deuteron 
break-up reaction. Notation as in fig. 2. 
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Fig. 4: The polarization transfer k for the exclusive proton - deuteron break-up reaction. 
Notation as in fig. 2. 
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